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The charge carrier mobility is a key parameter for the organic bulk heterojunction solar cell 
efficiency ft was recently shown that the interplay charge carrier transport and recombination, both 
depending on electron and hole mobilities, leads to a point of maximum power conversion efficiency 
at a finite mobility. Changes of bulk and surface recombination rate, however, can strongly influence 
this behavior. These processes were previously not considered adequately, as surface recombination 
velocities of infinity were implicitly assumed or bulk recombination parameters not discussed in 
detail. In this manuscript, using a macroscopic effective medium simulation, we consider how 
a reduced bulk recombination process in combination with finite surface recombination velocities 
affect the power conversion efficiency. Instead of a maximum efficiency at a specific charge carrier 
mobility, we show that with realistic assumptions and passivated surfaces the efficiency is increased 
further, saturating only at higher mobilities. Thus, a mobility optimisation is more important for 
the solar cell performance then previously shown. 

PACS numbers: 71.23.An, 72.20. Jv, 72.80.Le, 73.50.Pz, 73.63.Bd 



I. INTRODUCTION 

High power conversion efficiencies of organic solar 
cells are one key criterium for a use of organic semi- 
conductors in industrial applications. Until now, by 
improving process parameters, bulk heterojunction solar 
cell efficiencies of up to 5% are reached by several groups 
using well known materials like poly(3-hexylthiophcnc) 
(P3HT) and phen yl-C61-butyric acid methyl ester 
(PCBM) blendsPHJ Recently, by introducing new 
organic semiconductors providing advantageous energy 
levels and charge carrier mobilities, verified efficiencies 
of up to 6% and internal quantum yields of almost 100% 
could be achieved^. These examples show how impor- 
tant it is to design optimized organic semiconductor for 
solar cell applications. In order to predict the influence 
of charge carrier mobilities, numerical calculations have 
been performed, showing power conversion maxi ma at a 
certain mobility equal for electrons and hole d l 12 lFurthcr 
macroscopic simulations based on the detailed balance 
framework confirmed this prediction, but pointed out 
that a finite surface recombination at the electrodes can 
avoid a lowering of the power conversion efficiency at 
high charge carrier mobilities 9 . As simplification, all 
those calculation either implicitly assumed an infinite 
surface recombination rate for the contacts or do not 
discuss the influence of such details like the combined 
effects of finite surface recombination rates and a limited 
bulk, e.g. bimolecular, recombination. 

In this paper, we present how charge carrier mobili- 
ties and surface recombination processes of electrons and 
holes influence the power conversion efficiencies of an or- 



ganic solar cell. We do this by performing numerical 
simulations, also taking a physically reasonable reduced 
and capped bimolecular Langevin-like recombination into 
accounlP. 



II. SIMULATION 

To calculate current-voltage characteristics of or- 
ganic solar cells, a differential equation system of 
discretized Poisson and continuity equations for elec- 
trons and holes has to be solved as already shown 

For the active material the effective medium approach 
is used, based on a blend of organic semiconducting 
donor (polymer) and acceptor (fullerene) molecules 
under consideration of reachable donor acceptor in- 
terfaces within the diffusion length of charge carriers 
everywhere inside the active region. Due to the en- 
ergetic advantage excitions created e.g. by incident 
light on one molecule tend to separate to a coulom- 
bicly bound charge transfer state on a femto second 
timescale^. This leads to the model of an effective 
electrical band gap E g between the acceptor LUMO 
(lowest unoccumpied molecular orbital) and the donor 
HOMO (highest occupied molecular orbital) level, being 
able to describe the electrical properties of the whole 
device in terms an ambipolar effective medium with a 
band- like transport. For the numerical description, also 
input parameters such as the thickness L, the effective 
charge carrier densities N c , N v for HOMO and LUMO 
bands, the charge carrier mobilities for electrons fi n 
and holes /i p and the relative dielectric permeability e r 
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are required. Also the electron charge q, the vacuum 
permeability eo and the Boltzmann constant ks are used. 

On the active layer a predominantly electron con- 
ducting cathode and a predominantly hole conducting 
anode are added. Discontinuities of electrode work 
functions to the corresponding active material bands are 
considered by injection barriers <fr a , Furthermore, 
both contacts were assumed to possess finite surface 
recombination velocities. DDue to these contacts, the 
active material next to the cathode (anode) possesses 
a much higher electron (hole) charge carrier density, 
therefore, called the majority and a lower minority hole 
(electron) density. Consequently, within this paper, we 
denote the surface recombination, electrons crossing 
the cathode and holes crossing the anode interface, 
as the majority surface recombination Stdaji an d the 
vice versa, holes crossing the cathode and electrons 
crossing the anode interface, to be theminority surface 
recombination Smw- Reduced (finite) majority surface 
recombination velocities describe charges taking time to 
overcome injection or extraction barriers at the interface. 
For minority charges, a reduced surface recombination 
velocity results in a lower current of minority charges 
diffusing into the electrode, e.g. caused by mirror charge 
effects at the metallic electrode. In general real interfaces 
do not possess infinite recombination velocities. At the 
electrode minority charges recombine instantaneously 
with the present majority charges. Thus, a reduced 
minority surface recombination passivates the surface. 

At a given temperature T, we are now able to calculate 
the potential energy levels of Elumo, Ejjomo an d 
the corresponding quasi-Fermi levels of electrons E Fn 
and holes Ep v as well as the charge carrier density 
distribution inside the solar cell for electrons n and 
holes p. Energy levels and charge carrier densities are 
linked by the Boltzmann distribution, e.g. for electrons 
n = A c cxp (-q(E LUM o ~ E Fn )/k b T). 

Illumination of a solar cell can be described by adding 
a spatially constant charge carrier excitation G. This is 
opposed by a bimolccular, non-geminate recombination 
rate R, proven in disordered systems to accord to the 
theory of Langevin, describing recombination of c harges 
residing on different phases of the effective mediuni 11 ! 16 !. 
Recently, it was found that the experimental determined 
bimolecular recombinati on rate i s reduced in compari- 
son to Langevins 

> theor} ' 3 ' 5 ' 8 ' 15 -!^. We could explain the 
temperature dependence of this reduction factor by con- 
sidering the spatial charge carrier distribution gradients 
inside a device^. Nevertheless, the origin of the remain- 
ing constant reduction factor of £ = 0.1 is unknown so 
far. With this factor, the bimolecular recombination rate 
reads 



parameter 


symbol value 


unit 


temperature 


T 




300 


K 


effective band gap 


E g 




1.05 


eV 


relative dielectric constant 


e r 




4.0 




active layer thickness 


L 




100 


nm 


effective density of states 


N e , 


N v 


1.0 x 10 26 


l/m : 


generation rate 


G 




9.0 x 10 27 


l/m : 


injection barriers 






0.1 


eV 


Langevin reduction factor 


c 




0.1 





TABLE I. Parameters used for simulation. 



with 

7 = — (M« + A*p) (2) 

and the intrinsic charge carrier concentration 

m = ^/N C N V cxp (-qE g /2k b T) . (3) 

However, the reduced recombination rate is proved to 
be valid in low mobility materials only (< lcm 2 /Vs)P^. 
In terms of Eqn. faster charges would have a 

higher, unlimited recombination probability. Even so, 
as this assumption only holds for low mobilities, at high 
mobilities the bimolecular recombination is mobility 
independent, leading to a constant Langevin prefactor 
7. This can be compared with the mobility independent 
recombination rates of high mobility inorganic systems 
such as the compound semiconductor Cu(In,Ga)Se2. 

A possible but crude approximation simply cuts off the 
Langevin prefactor 7 at a certain point, now defined as 

7 = { i^T- + 0p) if + &>) ^ lent ^ 

\lcrit if (Mb + Mp) > la-it- 

Since j C rit is a physically reasonable but not directly 
accessible parameter, we convert it to a critical mobility 
fJicrit — ^^^Icritj assuming identical values for electrons 
and holes. In the following, we denote the loss according 
to Eqn. Q as capped Langevin recombination. A 
general overview of the underlying processes can be 
found irP. 

The presented calculation parameters listed in Ta- 
ble |TJ used for all presented calculations, were acquired 
by fitting the current-voltage characteristics of a rep- 
resentative 5% P3HT:PCBM solar cell (not shown) un- 
der assumption of infinite surface recombination veloc- 
ities. The charge carrier mobilities where determined 
experimentally 1 . 



R = (l (np - nf) , 
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FIG. 1. (Color online) Power conversion efficiency of a bulk 
hetero junction solar cell in dependence of charge carrier mo- 
bility equal for electrons and holes, for different surface re- 
combination velocities (solid lines). Under the assumption 
of a capped Langevin recombination (vertical lines at 1CP 8 
and 10 _4 m 2 /Vs), a finite minority surface recombination en- 
hances the solar cell efficiency to a saturated level at higher 
mobilities (dashed lines). Finite majority surface recombi- 
nation velocities generally lead to a reduction of the power 
conversion efficiency. 



III. RESULTS 

A. Variation of charge carrier mobility 

The calculated solar cell power conversion efficiencies 
for varied charge carrier mobilities, equal for elec- 
trons and holes, are plotted in Figure [T] The power 
conversion efficiency has an absolute maximum at a 
certain mobility (here around 10 _7 m 2 /Vs). Whereas the 
short circuit current increases with higher mobilities, 
the rather low power converison efficiencies at high 
mobilities are due to a drastical decrease of the open 
circuit voltage as already discussed irP. Taking finite 
minority surface recombination velocities of 10 _4 m/s, 
similar to the model of Scott and Malliaras^, this 
maximum is shifted to higher efficiencies. Taking the 
majority surface recombination rate additionally into 
account, we obtain a reduced power conversion efficiency. 

Introducing capped Langevin recombination rates 
in accordance with Eqn. Q, /i cr .;t set to 10~ 8 and 
10~ 4 m 2 /Vs, the power conversion efficiencies start 
to further increase with mobilities above the critical 
mobility, leading to a constant and saturated efficiency 
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FIG. 2. (Color online) Charge carrier mobility dependence 
of (top) open circuit voltage, (middle) short circuit current 
density and (bottom) fill factor, resulting in the power con- 
version efficiencies presented in Figure [l] The character- 
isics are shown for a reduced Langevin recombination (solid, 
( = 0.1) as well as for capped recombinations at mobilities 
of ficrit = 1CP 8 and 10~ 4 m 2 /Vs for different surface recombi- 
nation velocity combinations. In general, only the combina- 
tion of a finite minority surface recombination combined with 
a capped Langevin recombination rate shows an increase of 
Voc, Jsc and fill factor at mobilities higher than the critical 
mobility. 



at higher mobilities. Without finite minority surface 
recombination velocity, the power conversion efficiency 
is not influenced by the capped bimolecular recombi- 
nation. The reason becomes visible by considering the 
dependencies of V oc , J sc and FF on the charge carrier 
mobility (Fig. ^. 

The open circuit voltage shows a linear dependence 
on a logarithmic mobility axis, dropping from 0.8V at 
low mobilities (and low internal electrical fields) linearly 
to zero for high mobilities (and high internal electrical 
fields). Obviously, the difference between E g and V oc is 
the effective voltage drop due to the electric field inside 
the device. 

A combination of capped Langevin recombination — 
raising the internal charge carrier densities — and a 
finite minority charge carrier recombination — preventing 



charges to be extracted at the wrong contact — result 
in a higher and saturated V oc at mobilities above fi C rit- 
Without finite minority surface recombination, minority 
charges are extracted at both electrodes, leading to a 
reduced minority charge carrier density there. This 
effect has a strong effect on the band bending, enhancing 
the internal electric field. Thus, the open circuit voltage 
decreases. 

Beginning at low mobilities, the short circuit current 
rapidly increases to a saturated level for moderate 
and high mobilities. Under the influence of a finite 
minority surface recombination rate, this saturated 
current is further increased due to surface passivation, 
preventing minority charges to cross the electrode 
interface so that they finally recombine there with 
majority charge carriers. Introducing an additional 
majority surface recombination rate, creating a space 
charge at the electrodes, J sc is reduced again by en- 
hanced bulk recombination processes in those space 
charge areas. Furthermore, at high mobilities, s-shaped 
current-voltage characteristics reduce the short circuit 
current. A more detailed view on this effects will be 
discussed elsewhere^. Enabling capped bimolecular 
recombination rates, once again only the combination 
with an infinite minority surface recombination rate 
leads to a fast increase of short circuit current to a new 
saturated value at higher mobilities. Otherwise the short 
circuit current remains at a constant level or decreases. 

An analog behavior holds for the fill factor, showing 
a maximum around a mobility of 10 -7 to 10~ 4 m 2 /Vs. 
Exeeding the critical mobility suddenly increases the 
normally hill-shaped fill factor towards a saturation 
at higher mobilities, under the presence of a reduced 
minority surface recombination. 



B. Variation of surface recombination velocities 

At a charge carrier mobility of 10 _7 m 2 /Vs, the power 
conversion efficiencies are maximal (Figure [IJ dotted red 
line). For this constant mobility, the dependence of the 
power conversion efficiencies on recombination velocities 
is shown in Figure [3j 

Decreasing surface recombination velocities starting 
at infinity down to about 10 1 m/s do not influence the 
power conversion efficiency at all. In the following range 
down to 10 _2 m/s a dominant minority surface recombi- 
nation velocity results in an increase of power conversion 
efficiency until saturation is reached. For velocities 
lower than 10 _2 m/s, the majority surface recombination 
velocity starts to create space charges in the device 
leading to a drastic drop of power conversion efficiency. 
The behaviour of a capped Langevin recombination rate 
before and after the examined mobility of 10 _7 m 2 /Vs is 
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FIG. 3. (Color online) Intersection of Figure [TJ at a charge 
carrier mobility of 10 _7 m 2 /Vs, showing the influence of sur- 
face recombination velocities. A decreasing surface recombi- 
nation first improves the solar cell power conversion efficiency 
by a dominant minority surface recombination (if existent). 
At lower values, a possible majority surface recombination re- 
duces the efficiency drastrically. The behaviour for different 
Langevin capped mobilities remains identical. 

similar (dotted and solid lines). 

For a critical mobility of 10 _8 m 2 /Vs, V oc , J sc and 
fill factor are shown in dependence of the surface 
recombination velocity (Figure [4]). At the marked 
locations, the corresponding band structures for V oc and 
J sc conditions are shown in Figure [5] From these band 
structures the general influence of surface recombination 
velocities can nicely be seen as infinite values leading to 
an alignment of the quasi-Fermi levels and the electrodes 
work functions (Figure [4] a) . In contrast, finite velocities 
introduce a discontinuity of quasi-Fermi levels at the 
electrode interfaces as seen in (b) for minorities and (c) 
for majorities. The loss of minority charge carriers at 
interfaces can be identified by a high gradient of the 
quasi-Fermi levels, indicating a (recombination) current 
towards the electrode (J n oc n(x) ■ dEp n {x) / dx). 

The constantly improved power conversion efficiency 
by finite minority surface recombination velocities (Fig- 
ure [3]) can be traced back to a constant improvement of 
Voc, Jsc and fill factor (Figure |4| as soon as no essential 
minority recombination current is lost at the contacts 
(here around Smih = 10 _1 m/s). With a finite minority 
surface recombination, but without any perturbation by 
majority recombination, the device possesses maximal 
performance. This is shown in Figure pb). In contrast, 
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FIG. 4. (Color online) Surface recombination dependence of 
(top) open circuit voltage, (middle) short circuit current den- 
sity and (bottom) fill factor at different surface recombination 
velocities at a capped mobility of 10 _8 m 2 /Vs. With decreas- 
ing majority surface recombination, V oc as well as J 3C de- 
crease. In contrast, a finite minority surface recombination 
stabilizes V oc , the short circuit current density being domi- 
nated by majority surface recombination. The fill factor is 
also dominated by low majority surface recombination veloc- 
ities. 



a minority surface recombination current reduces the 
power conversion efficiency of a device without passi- 
vated contacts ( Figure [5^i). 

Starting from SMaj = 10~ 3 m/s, lower finite majority 
recombination velocities degrade the solar cell efficiency. 
As can be seen in Figure |5jc), the quasi-Fermi levels of 
the majority charge carrieres at the interfaces possess 
a discontinuity to the work function. Accordingly, 
charge carriers pile up inside the device, which increases 
the probability to recombine. Therefore, V oc starts 
to decrease linearly on a logarithmic scale (Figure 4]). 
Also the fill factor is reduced. At velocities smaller 
than Smuj = 10 _7 m/s, the short circuit current drops 
drastically. 

Bringing both surface recombination effects together 
(Figure |4j green line) leads to a constant and enhanced 
open circuit potential which is dominated by Sum- The 



FIG. 5. (Color online) Device energy structures at V oc (left 
column) and J sc (right column) for conditions marked in Fig- 
ure [4] Solid lines indicate the semiconductor LUMO (top) 
and HOMO (bottom) bands with the corresponding dotted 
quasi-Fermi levels for electrons (middle up) and holes (mid- 
dle down). A noninfinite surface recombination velocity intro- 
duces a jump of quasi-Fermi levels from the electrodes work 
function to the semiconductor quasi-Fermi levels. 



short circuit current as well as the fill factor first ben- 
efit from the positive influences of the minority surface 
recombination, but then are decreased by a dominant 
majority surface recombination, being even more pro- 
nounced in case of the fill factor. 



IV. CONCLUSION 

A lot of research efforts have been put into investi- 
gating the bulk recombination in bulk heteroj unction so- 
lar cells. Up to now, the effect of surface recombination 
has been largely ignored. Nevertheless, both types of re- 
combination need to be considered as they are crucial 
for optimizing the power conversion efficiency of organic 
photovoltaic devices. Previously, it has been shown that 
organic solar cells have an efficiency peak atmoderate- 
mobilities, but the efficiency decreases again, driven by 
a lowered open circuit voltage, at higher mobilities. In 
this contribution, we have shown that this behavior can 
only be observed for perfectly conducting contacts, im- 
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plying ideal charge extraction and Langevin recombina- 
tion, in which the decay rate is limited by the probability 
of charges finding one another for recombination. How- 
ever, both assumptions are not generally applicable. In 
organic devices, the charge extraction is not ideal, and 
Langevin recombination is only reasonable for rather low 
mobilities; for fast charge carriers, the recombination rate 
limiting step is not the finding, but the actual recombi- 
nation process. Therefore, in our calculations of power 
conversion efficiencies for bulk heteroj unction solar cells, 
we considered a bulk recombination following Langevin 
behaviour at low mobilities and a constant rate beyond. 
Furthermore, we accounted for surface passivating finite 
minority surface recombination velocities. The combi- 
nation of both physically justified processes leads to en- 
hanced power conversion efficiencies which do not de- 
crease once a maximum for a certain mobility has been 



reached, but increase further and finally saturate at high 
mobilities. Our results thus show that a mobility opti- 
misation can yield more in terms of performance then 
previously anticipated. 
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